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The first well founded perturbation theory for classical 
solid systems is presented. Theoretical approaches to ther- 
modynamic and structural properties of the hard-sphere solid 
provide us with the reference system. The traditional difficul- 
ties of all previous approaches are overcome. The perturba- 
tion is a first order term in an expansion of the Helmholtz free 
energy functional instead of an additive ad hoc term and the 
proper solid reference structure is used instead of some kind 
of mapping into the fluid structure at some effective density. 
As the theory reduces to the usual liquid perturbation theory 
in the uniform limit, it can describe consistently the complete 
phase diagram. Excellent results are obtained when applied 
to different model systems. 

PACS numbers: 64.60.-i, 64.70.Dv, 64.70. Kb 



In the last two decades, we have witnessed a continuous 
progress in the theory of non uniform classical fluids [Q,^ . 
Inside this context, where classical solids are regarded as 
an extremely non uniform system, a considerable effort 
has been devoted to develop a theory for solids. From 
this point of view, the theory would provide not only a 
description of solids but a unified view of classical sys- 
tems However, due to the inherent difficulties in the 
non uniform systems, this progress is still not comparable 
to that reached in uniform liquids. 

An important part of the enormous progress of the 
theory of simple uniform liquids has been due to the de- 
velopment of perturbation theories. Within this scheme, 
the liquid properties are described by those of a reference 
system (usually built with the repulsive part of the inter- 
action potential) modified by a perturbation (usually due 
to the attractive part of the interaction potential). The 
usefulness of these theories is based on the knowledge 
of the thermodynamic and structural properties of the 
reference system. The foremost, if not the unique, ref- 
erence system is the hard-sphere fluid, which describes 
many of the essential features of realistic interacting re- 
pulsive potentials. For this reason, much effort has been 
devoted to understanding its properties. Nowadays, the 
virial expansion, the analytically solvable Perkus-Yevick 
approximation and semi-empirical approximations sup- 
ply quite accurate results for thermodynamics and struc- 
ture up to the crystallization density |]l|. The seminal 
article by Andersen et al. Q establishes the relation be- 
tween the properties of a hard-sphere fluid and those of 
realistic repulsive interactions, completing the scheme of 



the perturbation theory for simple liquids. 

The same strategy has been pursued in the solid phase 
of classical systems. The aims of the theory are essen- 
tially two: the formal development of the theory itself 
and the determination of accurate thermodynamic and 
structural properties of the reference system. To this 
end, the density functional formalism has been success- 
fully elaborated to describe the thermodynamic proper- 
ties of the hard-sphere solid [Q. At present, there are 
different functional approaches for the Helmholtz free en- 
ergy which describe the solid phase and reduce to that of 
the fluid in the homogeneous limit Recently, it has 
been showed that the equation of state obtained from 
these functionals reproduces the simulation results quite 
well even up to almost close packing 1^. For all that, 
the thermodynamics of the hard-sphere system can be 
accurately described from low densities to almost close 
packing, including the crystallization phase transition. 

However, up to very recently, there was neither an ac- 
cessible theoretical approach to the structural properties 
of the hard-sphere solid nor a proper perturbation the- 
ory lH' -'^11 the perturbation approaches proposed by 
different authors, two of us included, are rather primi- 
tive 0-^. Basically, they follow the same scheme: the 
Helmholtz free energy of the solid is written as a sum 
of two terms, namely the reference and perturbation free 
energies. The former is assumed to be any of the available 
density functionals for the free energy of the hard-sphere 
solid while the perturbation is an ad hoc term built with- 
out any connection with the hard-sphere potential which 
is being assumed as the reference potential. Furthermore, 
the structure of the reference hard-sphere solid is substi- 
tuted in the perturbative term by some kind of mapping 
into the structure of the hard-sphere fluid at some effec- 
tive density. Typically, the Andersen et al. (WCA) cri- 
terium is used to divide the interacting potential, whereas 
the Barker-Henderson criterium ||l|,|l^ is used to fix the 
hard-sphere diameter of the reference solid. To worsen 
things, the Barker-Henderson criterium does not distin- 
guish uniform from non uniform systems. Any intent to 
do a consistent approach gives disastrous results. 

In spite of the crudeness and the fragility of these 
approaches, all pieces assembled achieved to describe 
the phase diagram of the Lennard- Jones system reason- 
ably well This success is not completely understood 
though some clues have been suggested Q . It seems that 
for large-ranged attractive interactions some kind of nu- 
merical error cancelation should occur. When they are 
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applied to systems with slightly more sophisticated po- 
tentials, they merely give a qualitative description of the 
thermodynamic properties or even fail completely. This 
is what happens when short-ranged attractive interac- 
tions are present. In any case, these theories are far from 
being satisfactory. 

The first step in order to develop a proper theory of 
perturbations for solids has been carried out recently by 
determining theoretically the average g{r) of the pair dis- 
tribution function p(^^(ri,r2) 1^. This average is defined 

by 
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(1) 



where V is the volume, p is the mean density, and dfl 
the differential solid angle aperture around ri2. Excellent 
results are obtained for the hard-sphere solid up to almost 
close packing. 

In this Letter, we develop a perturbation theory where, 
to first order, all the structural information needed is pre- 
cisely the function g{r). As it is usually done in pertur- 
bation theories for uniform liquids, we divide the inter- 
acting potential into the reference and the perturbative 
parts: (/){r) = (/)r{r) + 4>pif)- Then, we start from the 
general expression for the Helmholtz free energy of a non 
uniform system as a functional of the density p(r) which 
is written exactly as 0,|ll| 



(2) 



F[p{r)]^Fr[p{r)]+Fp[pir)], 



being Fr[p{r)] the Helmholtz free energy of the reference 
system at density p{r) and 

Fp[p{^)]^\J^ da j dvidY2P^'^\Yi,T2\a)(t>p{ri2), (3) 

where a is the coupling parameter for the interaction po- 
tential 4>(r]a) = (t)r{r) -)-a(/)p(r) and p^^^ (ri, r2; a) is the 
pair distribution function when the potential is 4'{r 12 ; a) 
but the density is p(r). To lowest order, Eq.(||) yields 

FpW)] = dridr2p(2)(ri,r2)</)p(ri2), (4) 

which after an appropriate integration becomes exactly 



Fp[p{v)] = 2TipN / drr-'gr{r)cj,p{r) 



(5) 



where N is the number of particles and g{r) is precisely 
the average of the pair distribution function given by 
Eq.(|l]). For any realistic system, the reference interact- 
ing potential is chosen to describe the rapidly varying 
short-ranged repulsive part of the interacting potential. 
However, the thermodynamic and structural properties 
of these systems are unknown. Therefore, a treatment to 
relate these properties to those of a hard-sphere solid is 



evidently necessary. The procedure is analogous to that 
of the theory of hquids. Lete(ri,r2) — exp[—P<f)\ri~T2\) 
be the Boltzmann factor. It depends only on ri2, but for 
clearness in some of the expressions below it is convenient 
to keep the formal dependence on the two locations ri 
and Y2. As in the liquid theory, let Ae(ri, r2) be the blip 
function, i.e., the difference of the Boltzmann factors of 
the reference potential and the hard-sphere potential of 
diameter dns, e^dri — r2|) — e^sdri — 1"2|)- If this diam- 
eter is comparable to the range of the reference potential, 
the blip function is different from zero in a small range of 
the order of S,dHs with ^ < 1. We can then expand the 
Helmholtz free energy functional of the reference system 
in powers of the blip function 

FrW)] = FhsW)] + 

- dvi dv2- -Ae ri,r2 

27 J de/fs(ri,r2) 



-higher-order terms. 



(6) 



The functional derivative of the Helmholtz free energy 
with respect to the Boltzmann factor is easily obtained 



(7) 



where the y(ri,r2) function is the generalization of the 
y{ri2) function in uniform systems. This last one is de- 
fined by e~'^'^^''i^^(?(ri2), whereas the former is defined by 

■ Introducing Eq.(0) into Eq.(|) and 
after some simple algebraic steps, it is found that 

FrW)]^FHs[p{r)] + 

27: pN j (irr^yHs(?')Ae(r) -|- higher-order terms, (8) 

where ynsif) is the average of the y^rs (ri,r2) func- 
tion which can also be expressed as yHs{f) = 
exp{P(l)Hsif))gHs{i-)- Until here, the hard-sphere diam- 
eter has remained arbitrary. Now we specify it by impos- 
ing that the first order term in the functional expansion 
Eq.(p|) be zero. This yields 



drr yHs{r) = / drr yHs{r)exp{(34>r{r)). (9) 



Notice that it is analogous to the well known WCA cri- 
terium |^ in liquid theory except that the y{r) function 
instead of the y{r) function is used. From Eq.(^ and 
Eq.(^ and after averaging, we find 



yr{T) = VHsir) + higher-order terms. 



(10) 



which can be used to lowest order to evaluate the pertur- 
bation term Eq.(||). It is interesting and easy to prove 
that, if Eq.(9) holds, the convergence properties of the 
expansions (3) and ( p^ are similar to those found in the 
uniform limit ||l^ , i.e. 
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2npN J drr^yHsir)Aeir) + 0(e*), (11) 

and 

yrir)=yHsir) + Oie)- (12) 

With the above expressions the theoretical formahsm 
is completed. It is straightforward to see that in the 
uniform limit the theory reduces to the WCA perturba- 
tion theory of simple liquids. Observe also that all the 
problems of previous approaches are automatically over- 
come: the formalism provides a proper expansion of the 
free energy with precise thermodynamic and structural 
properties of a reference solid, namely any of the func- 
tionals for the Helmholtz free energy and the correlation 
function g{r), and the criterium to determine the appro- 
priate reference system enhances the convergence of the 
expansion. 

We have applied the theory to the Lennard- Jones (LJ) 
system and to an extremely short-ranged square well po- 
tential. In Table I we show the LJ liquid {pi) and solid 
(ps) densities at coexistence at several temperatures (all 
in LJ parameters units) obtained from the present the- 
ory, using Tarazona functional ||l^ for the reference hard- 
sphere solid, compared to simulation data of Hansen and 
Verlet |lj]. The theoretical predictions are quite good 
and the small deviations are quantitatively similar to 
those obtained by previous ad hoc approaches. Never- 
theless, contrary to what it is predicted by the later, the 
present theory gives the correct slopes of the coexistence 
densities as functions of the temperature. This improve- 
ment is undoubtfuUy due to the adequate dependence of 
the hard-sphere diameter with the density. Table I also 
shows the solid Lindemann parameter at coexistence. It 
is significantly smaller than the simulation result. This is 
due to the functional approach we have used for the refer- 
ence hard-sphere solid, as it is well known. If any other 
functional were used, the Lindemann parameter would 
improve but no significant changes would be detected for 
the energies and, therefore, for the coexistence densities 
shown in Table I. However, the important point is that, 
for the first time, a classical system can be studied with 
a unique, proper and consistent theory. 

Much more impressive are the results obtained by the 
theory when it is applied to a system with an extremely 
short-ranged attractive interaction. Simulation results 
have recently proved that this kind of systems present 
quite interesting isostructural solid-solid transitions. Up 
to now, as one could expect, no theory had been able 
to give a reasonable quantitative account of this behav- 
ior |p^ , |l7| . Therefore, these systems offer an interesting 
and stringent test for the present theory. Figure 1 shows 
the phase diagram of the square well potential charac- 
terized by a short width S/dns ~ 0.02. The dashed 



lines of this figure, which corresponds to simulation re- 
sults jl^, show the coexistence densities of the liquid- 
solid transition and, inside of the solid region, the coex- 
istence curve of the solid-solid transition. Observe that, 
due to the weakness of the attractive potential, there is 
no fluid condensation. However, the solid-solid transi- 
tion is a kind of solid condensation which arises from a 
different mechanism than the usual fluid condensation. 
It is related with the commensurability of the lattice pa- 
rameter of the solid structure with the well width of the 
potential [ p^Jl6| . The dotted lines correspond to the 
perturbation weighted density approximation (PWDA) 
P,p^ which was, to our knowledge, the best approach so 
far able to describe the complete phase diagram of this 
system. The continuous lines correspond to the present 
theory using Tarazona functional [|l^ for the reference 
hard-sphere solid. The results using any other functional 
would be indistinguishable in the scale of this figure. The 
dramatic improvement is clear and the global agreement 
of the theoretical predictions with the simulation results 
is quite good. Equivalent quantitative agreement is ob- 
tained for the phase diagrams corresponding to square 
wells of different widths. 

In summary, we have developed a perturbation the- 
ory with remarkable characteristics. It is the first well 
founded perturbation theory for classical solids. It re- 
duces to the well known WCA perturbation liquid theory 
in the appropriate uniform limit. Its versatility allows its 
use with any of the functional approaches available for 
the Helmholtz free energy of hard spheres. It gives quite 
good results even for systems where other ad hoc ap- 
proaches have failed. Finally, but not less important, the 
computational effort is drastically reduced. The cum- 
bersome mapping of the structure of the solid into the 
structure of the liquid at certain effective density, which 
usually involves tedious autoconsistent and recursive pro- 
cesses, is replaced by a simple integral. 
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FIG. 1. Phase diagram of the square well system for the 
well width S/dns = 0.02. The temperature is given in well 
depth units. The solid line is the prediction of the present 
theory. Dashed line corresponds to the simulations by Bolhuis 
and Frenkel. Dotted line is the prediction of the PWDA. 

TABLE L Lennard- Jones liquid {pi) and solid (ps) 
densities at coexistence at several temperatures (all in 
Lennard-Jones parameters units) as predicted by the simu- 
lations of Hansen and Verlet and the present theory. The 
Lindemann parameter L of the solid phase at coexistence is 
also shown. 
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